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Abstract 

We propose a simple approach to the calculation of the entanglement entropy of a spherically 
symmetric quantum system composed of two separate regions. We consider bound states of the 
system described by a wave function that is scale invariant and vanishes exponentially at infinity. 
Our result is in accordance with the holographic bound on entropy and shows that entanglement 
entropy scales with the area of the boundary. 
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I. INTRODUCTION 



Entanglement is one of the fundamental features of quantum mechanics and has led to the 
development of new areas of research, such as quantum information and quantum computing. 
Some reviews and recent results on entanglement entropy in many-body systems, conformal 
field theory and black hole physics can be found in 

Einstein, Podolsky and Rosen (EPR) proposed a thought experiment 6| to prove that quan- 
tum mechanics predicts the existence of "spooky" nonlocal correlations between spatially 
separate parts of a quantum system, a phenomenon that Schrodinger [?| called entanglement. 



Afterward, Bell [8j] derived some inequalities that can be violated in quantum mechanics but 
must be satisfied by any local hidden variable model. It was Aspect (9] who first verified 



in laboratory that the EPR experiment, in the version proposed by Bohm [lOj, violates 
Bell inequalities, showing therefore that quantum entanglement and nonlocality are correct 
predictions of quantum mechanics. 

A renewed interest in entanglement came from black hole physics: as suggested in ll|, [12] , 
black hole entropy can be interpreted in terms of quantum entanglement, since the horizon 
of a black hole divides spacetime into two subsystems, such that observers outside cannot 
communicate the results of their measurements to observers inside, and vice versa. Black 
hole entanglement entropy turns out to scale with the area A of^ the event horizon, in 
accordance with the renowned Bekenstein-Hawking formula |l3 . 
where Sbh is the black hole entropy and Ep is the Planck length. 
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Let us consider a spherically symmetric quantum system, composed of two regions A 



T]). The entanglement entropy 
for many-body systems. This 



and B separated by a spherical surface of radius R (Fi, 
of each part obeys an "area law", as discussed e.g. in [2 
result can be justified by means of a simple argument proposed by Srednicki in 12|. If we 
trace over the field degrees of freedom located in region B, the resulting density matrix p A 
depends only on the degrees of freedom inside the sphere, and the associated von Neumann 
entropy is Sa = — Tr(p A lnp A ). If then we trace over the degrees of freedom in region A, we 
obtain an entropy Sb which depends only on the degrees of freedom outside the sphere. It is 
straightforward to show that Sa = Sb = S, therefore the entropy S should depend only on 
properties shared by the two regions inside and outside the sphere. The only feature they 
have in common is their boundary, so it is reasonable to expect that S depends only on the 
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FIG. 1: A quantum system composed of two parts A and B, separated by a spherical surface of 
radius R. 



area of the boundary, A = AixR 2 . 

The area scaling of entanglement entropy has been investigated much more in the context 
of quantum field theory than in quantum mechanics. In order to bridge the gap, in this 
paper we study the entanglement entropy of a quantum system composed of two separate 
parts (Fig. [1]), described by a wave function ip, which we assume invariant under scale 
transformations and vanishing exponentially at infinity. In Section II I II we will show that 
the entropy S of both parts of our system is bounded by S < 77^, where 77 is a numerical 
constant related to the dimensionless parameter A appearing in the wave function ip. This 
result, obtained at the leading order in A, is in accordance with the so-called holographic 



13, [l8|, ll9|: S < where A is the area of 



bound on the entropy S of an arbitrary system 
any surface enclosing the system. 

In Section [III we present the main features of our approach, focusing in particular on 
the properties of entanglement entropy and on the form of the wave function describing the 
system. In Section II I II we calculate analytically the bound on entanglement entropy. In 
Section IIVI we summarize both the limits and the goals of our approach. 

II. MAIN FEATURES OF THE MODEL 

Let us suppose that a quantum system consists of two parts, A and B, which have 
previously been in contact but are no longer interacting. The variables g, r describing the 
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system are subjected to the following constraints (see Fig. [I]): 

< g < R region A 
r > R region B , 

where R is the radius of the spherical surface separating the two regions. It is convenient to 
introduce two dimensionless variables 

X= R> y= R> (1) 
subjected to the constraints < x < 1 and y > 1. 

In the following we will assume that the system is spherically symmetric, in order to treat the 
problem as one-dimensional in each region A and B, with all physical properties depending 
only on the radial distance from the origin. 



A. Entanglement entropy 

Let ip(x, y) be a generic wave function describing the system in Fig. (TJ composed of two 
parts A and B. 



As discussed e.g. in 
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211 ]. we can provide a description of all mesauraments in region A 



through the so-called density matrix p A (x, x')\ 

p A (x, x') = / d 3 y?p*(x, y)tp{x', y) , (2) 



where d 3 y is related to the volume element d 3 r in B through the relation d 3 r = R 3 d 3 y, with 
d 3 y = y 2 sin 9d9 d<p dy in spherical coordinates. 

Similarly, experiments performed in B are described by the density matrix p B (y, y'): 

„,(*,*')= vW*rS), P) 



where d 3 x is related to the volume element d 3 g in A through the relation d 3 g = R 3 d 3 x, with 
d 3 x = x 2 sin 6d9 d(f) dx in spherical coordinates. 

Notice that p A is calculated tracing out the exterior variables y, whereas p B is evaluated 

tracing out the interior variables x. 

Density matrices have the following properties: 

1. Trp = 1 (total probability equal to 1) 



2. p = p^ (hermiticity) 

3. pj > (all eigenvalues are positive or zero). 

When only one eigenvalue of p is nonzero, the nonvanishing eigenvalue is equal to 1 by virtue 
of the trace condition on p. This case occurs only if the wave function can be factorized into 
an uncorrelated product 

V>(ar, y) =ip A (x) -ip B (y)- (4) 

States that admit a wave function are called "pure" states, as distinct from "mixed" 
states, which are described by a density matrix. A quantitative measure of the departure 
from a pure state is provided by the so-called von Neumann entropy 

S = -Tr(plogp) . (5) 

5* is zero if and only if the wave function is an uncorrelated product. 

The von Neumann entropy is a measure of the degree of entanglement between the two parts 
A and B of the system, therefore it is called entanglement entropy. 

When the two subsystems A and B are each the complement of the other, entanglement 
entropy can be calculated by tracing out the variables associated to region A or equivalently 
to region B, since it turns out to be S A = S B . 

B. Wave function 

As already said, the spherical region A in Fig. [T]is part of a larger closed system A U B, 
described by a wave function ip(x, y), where x denotes the set of coordinates in A and y the 
coordinates in B. 

We will exploit for ip the following analytic form: 

^y) = Ce~^\ (6) 

where C is the normalization constant and A is a dimensionless parameter, that we assume 
much greater than unity (A ^> 1). 

If the system is in a bound state due to a central potential, the complete wave function ip 
should contain an angular part expressed by spherical harmonics Yi m (8, <p): 

4>(x, y; 9, 4>) = C Y lm (6, <p) ip(x, y) , 
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where C is the normalization constant. If the angular momentum is zero, the angular 
component of the wave function reduces to a constant Yqq(8, (f>) = which can be 

included in the overall normalization constant C = C / \^Att appearing in the expression (J6]) 
of the wave function if). 

In order to justify the form (J6j) of the wave function if), let us list the main properties it 
satisfies. 

1) ip depends on both sets of variables x, y defined in the two separate regions A and B, but 
it is not factorizable in the product (j4j) of two distinct parts depending only on one variable: 

ip(x, y) ^ if) A (x) ■ ip B (y) . 

This assumption guarantees that the entanglement entropy of the system is not identically 
zero. 

2) if) depends on the variables x, y through their ratio y/x, hence it is invariant under scale 
transformations: 

x — > fix and y — > fiy, with /i constant . 

3) if) has the asymptotic form of an exponential decay. 

The last ansatz on if) is equivalent to consider the quantum state of a central potential 
vanishing at infinity, with negative energy eigenvalues. The asymptotic form of the radial 
part /(r) of the wave function describing this state is 

f{r) = C B e-" = C B e~ KR y, (7) 

where r —>■ oo is the radial distance from the origin, C B is the normalization constant and 
y = r/R. In particular, we could consider a particle with mass m and negative energy E, in 
a bound state due to a central potential going to zero as r — > oo. In this case, the asymptotic 
form of the Schrodinger equation, in spherical coordinates, is given by 

d 2 f(r) o r/ N , j2m\E\ 

-jy = K 2 f(r), with «=V_J-i. (8) 

Apart from the normalization constant, f(r) coincides with the restriction if> B {y) of the wave 
function if)(x, y) to the exterior region B, as seen by an inner observer localized for instance 
at x = 1, i.e. on the boundary between the two regions: 

1>M = il>(x,v)\ KSsl = Ce-». (9) 
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By comparing the asymptotic behaviour of the wave functions © and as y — > oo, we 
find: 

R 

A = 7 — , with j = k£p, (10) 
ip 

where the Planck length £p = (hG/c 3 ) 1 ^ 2 has been introduced to make the parameter 7 
dimensionless, without any further physical meaning in this context. In Section HTT1 we will 
assume A ^> 1, which is always true in a system with R sufficiently larger than £p, as it 
easily follows from Eq. (flQj) . 

If the inner observer is not localized on the boundary but in a fixed point xq inside the 
spherical region (with < x < 1), then the expression (TTU]) of A has to be multiplied by x . 
Notice that the dependence of A on the radius R of the boundary has been derived by 
imposing an asymptotic form on the wave function ip(x, y) as y — ► 00, with respect to a 
fixed point x ~ 1 inside the spherical region of the system in Fig. [IJ 

In Section II III we will show that the entropy of both parts of our system depends on 
A 2 ~ R 2 /£p, i.e. on the area of the spherical boundary. The area scaling of the entanglement 
entropy is, essentially, a consequence of the nonlocality of the wavefunction if)(x, y), which 
establishes a correlation between points inside (x ~ 1) and outside (y — ► 00) the boundary. 



III. ANALYTIC RESULTS 



We normalize the wave function ([6]) by means of the condition 

d 3 x / d 3 yi>*(x, y)i/)(x, y) = 1 , 



with d 3 x = x 2 sin 9d6 d<fi dx and d 3 y = y 2 sin QdQ dip dy in spherical coordinates. Under the 
assumption A> 1, the normalization constant C turns out to be 

C*2\^. (11) 

Let us focus on the interior region A represented in Fig. [U We calculate the density matrix 
p A (x, x 1 ) by tracing out the variables y related to the subsystem B, as expressed in Eq. (j2j): 

p A (x, x') = / d 3 yip*(x, y)ip(x', y) 

J B 

C 2 xx \ x+x' 

w 4vr- e ~ x ^ r , (12) 

A x + x' v ' 
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where we have assumed A > 1. 

It is easy to verify that the density matrix p A (x, x') satisfies all properties listed in Section 

EH 

1. Total probability equal to 1: 

d 3 xp A (x, x) = 1 Tr(p A ) = 1 . 



2. Hermiticity: 

p A (x, x') = p* A (x', x) p A = p\. 

3. All eigenvalues are positive or zero: 

P A {x, x')>o y x ,x'e(o,i)^(p A ) J >o 

The entanglement entropy (J5]) can be expressed in the form 

J A 

Substituting the expression f[T2"l) of p A (x, x'), with x' = x, we find: 
S A « 4vr— 

We can maximize the previous integral by means of the condition 

£ -2X/x < g -2A V X G (0, 1) . 

The entanglement entropy Sa turns out to be bounded by 



f d 3 xe~ 2X/x 


A log 


(a ° 2 


x\ 


' A 


2 & 


ri 


V. 



Sa < (4tt) 2 C 2 e~^{\ - 1 [a log ( 4 vr^) + ft] } 



C 2 



(13) 



with a = | and b = -^(1 + 4 log 2). 

By inserting the expression ( TTTT) of the normalization constant C, we obtain 

1 f 1 2 r 
Sa< -A 2 |l- y [alog(A/7r) + 6 

If we neglect the sub leading terms in A ^> 1 and substitute A = jR/ip, as given in Eq. 
(flOjh we finally find: 

^4 1 

5a ~^4£2"' with V = ^-l , ^ 
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where A = AttR 2 is the area of the spherical boundary in Fig. [TJ The result (JT 



accordance with the holographic bound on entropy S < ^4-, introduced in 17], |l8], and 



is m 



shows that the entanglement entropy of both parts of our composite system obeys an "area 



law", as discussed e.g. in 



2| for many-body systems. 
For a particle with energy E and mass m, satisfying the asymptotic form (JSJ) of the 
Schrodinger equation, we can express the parameter 77 in the form 

2 m\E\ , . 

Sir mpc z 

where we have combined Eqs. ([8]), (flQj) . ({T41) and have introduced, for dimensional reasons, 
the Planck mass mp = (hc/G) 1 ^ 2 . Under the assumptions \E\ <C mpc 2 and m < mp, we 
obtain ?j < 1, therefore in this case the bound (jl~4"|) on entropy is much stronger than the 
holographic bound S < 

All calculations performed in this Section can be repeated focusing on the exterior region 
B represented in Fig. CD By tracing out the interior variables x, as in Eq. ([3]), the density 
matrix p B (y, y') turns out to be 

P B (y^y') = / d 3 x?p*(x, y)ip(x, y') 

J A 

Q2 e -\{y+y') 

w 4vr- -, (16) 

A y + y' 

where we have substituted the expression of the wave function ip and have applied the 
usual assumption A 3> 1. 

Analogously to Eq. ( [131) . the entanglement entropy is given by the integral 

S B = - [ d 3 yp B (y, y)log[p B (y, y)\ . (17) 

J B 

The numerical evaluation of the previous integral confirms that the entropy bound calculated 
tracing out the interior variables x coincides, under the assumption A ^> 1, with the entropy 
bound evaluated tracing out the exterior variables y, i.e. Sa = Sp- 

IV. CONCLUSIONS 

We have proposed a simple approach to the calculation of the entanglement entropy of a 
spherically symmetric quantum system. The result obtained in Eq. (Tl4l . Sa ^ V -$r, is in 
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accordance with the holographic bound on entropy introduced in {17I . [Q and with the "area 
law" discussed e.g. in [J Q|. Our result, in fact, shows that the entanglement entropy of 
both parts of the system in Fig. [TJ depends only on the area of the boundary that separates 
the two regions. 

The area scaling of the entanglement entropy is a consequence of the nonlocality of the wave 
function, which relates the points inside the boundary with those outside. In particular, we 
have derived the area law for entropy by imposing an asymptotic behaviour on the wave 
function i[>(x, y) as y — > 00, with respect to a fixed point x ~ 1 inside the interior region. 
The main limit of our model is that we have considered only one particular form of the wave 
function ip. However, more general forms of ip might be considered, hopefully, in future 
developments of the model. Let us finally stress that our results are valid at the leading 
order in the dimensionless parameter A> 1 appearing in the wave function tp of the system. 

The treatment presented in this paper for the entanglement entropy of a quantum system 
is an extremely simplified model, but the accordance of our result with the holographic bound 
and with the area scaling of the entanglement entropy indicates that we have isolated the 
essential physics of the problem in spite of all simplifications. 
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